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Observing scattered X-ray radiation from gamma-ray bursts: a 
way to measure their collimation angles 
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Abstract. An abstract should be given 



here are observational facts and theoretical arguments 
for an origin of gamma-ray bursts (GRBs) in molecular 
clouds in distant galaxies. If this is true, one could detect 
a significant flux of GRB prompt and early afterglow X- 
ray radiation scattered into our line of sight by the molec- 
ular and atomic matter located within tens of parsecs of 
the GRB site long after the afterglow has faded away. The 
scattered flux directly measures the typical density of the 
GRB ambient medium. Furthemore, if the primary emis- 
sion is beamed, the scattered X-ray flux will be slowly 
decreasing for several months to years before falling off 
rapidly. Therefore, it should be possible to estimate the 
collimation angle of a burst from the light curve of its 
X-ray echo and a measured value of the line-of-sight ab- 
sorption column depth. It is shown that detection of such 
an echo is for the brightest GRBs just within the reach of 
the Chandra and XMM-Newton observatories. 
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1. Introduction 

One of the fundamental questions relating to the gamma- 
ray burst (GRB) phenomenon is how much energy is re- 
leased during a GRB. Given a measured burst fluence, 
the inferred energy release is proportional to Oq/2, where 
9q is the opening angle of the coUimated relativistic fire- 
baU (Piran 1999) that is probably smaU {0o < 1). Apart 
from constraining the burst energetics, 6q also determines 
the event rate of GRBs in the universe (which is propor- 
tional to Oq'^). Thus, the beaming angle is a very impor- 
tant quantity in GRB research. 

There have been attempts to estimate 6*0 via observ- 
ing a break in GRB afterglow light curves and interpreting 
this break as due to a slowing down of the jet to F < 



Send offprint requests to: sazonov@mpa-garching.mpg.de 



(here F is the bulk Lorentz factor of the jet). Summarizing 
the information on jet break times in more than a dozen 
of GRBs Frail et al. (2001) have inferred values for 
ranging from 1° to more than 25°. These authors came 
to the conclusion that the gamma-ray energy release is 
narrowly clustered around 5 x 10^° ergs. However, since 
their derivation was based on a particular interpretation of 
the observed breaks and that other explanations for such 
breaks exist (as mentioned by Frail et al. 2001), these mea- 
surements of 6*0 should necessarily be considered tentative 
for the time being. 

Another fundamental issue is the environment in which 
GRBs occur. High gas column densities (iVn ~ 10^^ cm~^) 
toward GRB locations have been inferred from a spec- 
tral analysis of a sample of X-ray afterglows observed 
with BeppoSAX (Owens et al. 1998; Galama et al. 2001; 
Reichart & Price 2002). Such values of Nn are typical of 
Galactic giant molecular clouds and therefore a strong in- 
dication that GRBs occur in star forming regions of their 
host galaxies, which in turn argues for collapse of massive 
stars as their sources (Woosley 1993; Paczynski 1998). 

In this paper we propose an observational test that 
enables direct determination (with some reservations) of 
both the collimation angle of a GRB and the typical den- 
sity of the medium surrounding it on scales of pes to tens 
of pes typical of molecular clouds and complexes. This 
method consists of observing the location of a bright GRB 
with a powerful X-ray telescope several times during the 
first months to years after the burst. We predict that in 
such observations a weak X-ray flux may be detected, 
which will be radiation from the burst reaching us, de- 
layed by scattering from the ambient medium. 

The idea at the basis of the proposed method - to 
search for scattered GRB radiation - is not new. The 
theme of scattered X-ray emission was first introduced by 
Dermer, Hurley & Hartmann (1991) in the context of the 
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then popular model of the galactic stellar binary origin of 
GRBs. Madau et al. (2000) considered Compton echoes 
from gamma-ray bursts arising in the circumburst envi- 
ronment within a fraction of a pc of the GRB site. These 
authors focused on the case where the primary burst emis- 
sion is coUimated away from us. Ramirez-Ruiz et al. (2001) 
furtherelaborated this scenario. Esin & Blandford (2000) 
studied the scattering of GRB early optical emission on 
the surrounding dust, while Meszaros & Gruzinov (2000) 
considered a similar phenomenon - small-angle scattering 
of X-rays by dust. 

GhiscUini et al. (1999) and Bottcher et al. (1999) have 
computed the time dependence of fluorescence emission 
in the iron Ka line resulting from the interaction of the 
GRB radiation with the surrounding matter under the as- 
sumption that the burst emission is isotropic. Detection 
of such a fluorescent line would make it possible to de- 
termine the redshifts of GRBs with invisible optical after- 
glows. However, the predicted flux in the line is apparently 
too low to be detected with the current generation of X- 
ray telescopes unless the density of interstellar material at 
the GRB site is very high {n ^ 10'* cm""^). We note here 
(see that the scattered X-ray flux in the relevant 0.3- 
5 keV band will typically be an order of magnitude higher 
than the flux of fluorescence emission in the same band. 
As we shall see, this enhancement factor proves crucial 
for the problem at hand, as it makes the emergent signal 
detectable with the present-day X-ray telescopes provided 
that GRBs indeed originate within molecular clouds. 




Fig. 1. Schematic diagram (cross section) of the model in 
which a coUimated GRB occurs at the center of a molec- 
ular cloud and the observer sees GRB X-ray radiation 
scattered by the gas. The scattering sites contributing to 
the detected flux at a given time lie on a constant-delay 
paraboloid between the jet and cloud boundaries. 



2. Delayed scattered X-ray emission from GRB 

We base our treatment on a simple model depicted in 
Fig. 1^. A GRB occurs at the center of a spherical cloud 
of gas of constant density n (equivalent number density of 
hydrogen atoms; the gas may be a mixture of atoms and 
molecules) and radius R. The host galaxy is located at 
a redshift z. The blastwave expands into a cone with an 
opening angle so that only observers located within this 
cone can directly receive X-ray and gamma-ray radiation 
from the burst. Note that since the bulk Lorentz factors of 
GRB fireballs are believed to be very large (F ^ 10^-10"^), 
observers in other directions will detect very low fluxes 
of radiation and at energies below X-rays. In our model, 
we are located exactly on the symmetry axis of the jet. 
Another parameter of the model is the GRB fluence, ^x, 
in some specified X-ray band £2]- 

The scattered X-ray emission is observed from Earth at 
a time t after the GRB. The scattering surface of constant 
delay will be a paraboloid with its focus at the GRB site 
and its axis along our line of sight, 
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(Blandford & Rees 1972), where is the scattering angle 
and c is the speed of light. The factor (1 + z)^^ takes into 
account cosmological time dilation. As is detailed below 
(in the fluence of the early X-ray afterglow [t <; lO'* s) 



may for some GRBs constitute a significant fraction of the 
X-ray fluence of the burst itself. The usual interpretation 
for such early afterglow emission is that it arises from the 
interaction of the jet with the external medium at a dis- 
tance of 0.01-0.1 pc from the site of the burst. Therefore, 
equation (|l|) will remain an appropriate approximation 
when considering the echo associated with the early after- 
glow if i? ^ 1 pc, which we assume to be the case. 

Only positions (r < R,0 < 60) will contribute to the 
scattered X-ray flux. The corresponding range of scatter- 
ing angles is < 9 < Oq, where the critical angle (see 
Fig. 0) 
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The scattered X-ray flux in [Ei, E2] is then 
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da dr 
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where dfi = 27rdcos 6 and dr/dt = c(l -I- z)~*(l — cos9)~^. 
The differential cross section for scattering, which in- 
cludes coherent (Rayleigh) and incoherent (Raman and 
Compton) scattering, is given by 



(4) 
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Fig. 2. Scattered X-ray flux as a function of elapsed time 
after a GRB at z = 1 with = lO'^ erg cm-^. The 
different sohd hues correspond to a GRB origin in the 
center of a uniform molecular cloud with i? = 10 pc and 
n = 10'^ cm~'^, and different burst coUimation angles: 
9o = TT, 25° and 15°. Also shown are the light curves for 
a GRB with 6*0 = 25° located at a; = 0.5i?, z = -0.5i? 
(see Fig. ^ in the same uniform cloud (short-dashed line) 
and in a cloud of radius R ^ 10 pc with a density law 
n{r) = no{R/r), uq — 5 x 10^ cm^"^ (long-dashed line). 
The horizontal dash-dotted line represents the scattered 
signal from the counterjet. The dotted line represents X- 
ray afterglow emission with parameters similar to the af- 
terglow of the bright GRB 000926 observed with Chandra 
at a late time t ~ 10^ s. 



where (Tt is the Thomson cross section. The coefficient 
A{9) allows for the fact that Rayleigh scattering of X- 
rays through small angles on molecules of hydrogen and 
atoms of heavier elements, most importantly helium, is 
more efficient (calculated per electron) than scattering on 
atomic hydrogen. A{9) takes values between 1 and 2 for 
— > 0, depending on the fraction of molecular hydrogen 
and helium in the interstellar medium, and is close to 1 for 
large scattering angles (e.g. Sunyaev & Churazov 1996). 
As we are primarily interested in situations where the jet 
opening angle is small (^ 30°) and the scattering medium 
is a molecular cloud, we shall adopt A{d) ^ 1.5 for our 
estimates below. 

On integrating over the solid angle in equation (^ we 

get 



Fx = -A5xn(l + z)-icaT/(^o,0i) 



f{0o,0i) = 



(cos do — COS 6i ) (cos 00 + cos 6'i + 2) 
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The function /(^Oj ^i) carries information on the time de- 
pendence of the scattered flux. Introducing a critical time 

R{l + z){l-coseo) 
ti = , (7j 



we may rewrite equation (p|) as 



/ = 21n 



1 — cos 6*0 
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2\nj fortcRB < i < ^i, 
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where icRB is the duration of the burst. Wc notice that 
the scattered flux tends to infinity as i — > 0, being lim- 
ited only by the finite GRB duration. This reflects the 
fact that at early times the dominant contribution to the 
flux comes from small scattering angles while the scat- 
tering front sweeps rapidly through the cloud toward us: 
dr/dt oc 9^^ cx t^^. We also see that if the GRB emission 
is coUimated in a narrow cone, then the slow (logarithmic) 
decline of Fx will be followed by a rapid drop to zero when 
the critical angle 9i{t) approaches the jet opening angle 
9o. The scattered flux will vanish completely at t = ti. 

Equations (||) , (0) and (||) provide a complete descrip- 
tion of the scattering effect. Substituting some values for 
the parameters, we get 



Fx = 



for 



t < ti 



5 X 10" 
A 



erg cm 
Sx 



1 fit) 



1.5 10 ^ erg cm ^lO'^cm 1 



10^- 



R 



■(1 -I- z)(l - cos 6*0), 



(9) 



(10) 



(5) 



10 pc^ 

and Fx = for t > ti. Note that / ~ 4 corresponds to 
t - O.lti. 

Fig. H shows examples of light curves of scattered GRB 
X-ray emission. One of the light curves represents the echo 
from a possible counterjet expanding away from us; the 
corresponding flux is very weak. Also plotted is the ex- 
trapolated light curve of the X-ray afterglow of the bright 
GRB 000926. This afterglow was observed with Chandra 
two weeks (i ~ 10^ s) after the GRB, when the flux was 
still decaying as a power-law with an estimated index 
a ^ —2 (Piro et al. 2001). This is the latest performed 
measurement of a GRB X-ray afterglow to our knowl- 
edge and as such it may be indicative of how steep the 
decay of X-ray afterglows is at t ;> 10^ s. Clearly, obser- 
vations aimed at detecting the scattered radiation should 
be scheduled for a not too early moment of time, so that 
the afterglow signal has faded away. It appears that wait- 
ing for a few months after the GRB should typically be 
enough. 
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Suppose now that it is possible in an experiment to 
measure ^x, z, Fx{t) and ti. One will then be able to 
make the following estimates: 

90 = arccos (l - 
and 

1 5 

n = 2 X lO^cm"^ — (1 + z) 

10-5 erg fx ft) 4 

~^ Sx 10-15 erg cm-2 s-i /W 

with /(i) being given by equation (||). 

It is evident that by performing one or two measure- 
ments of the scattered flux at well separated times t ^ti, 
i.e. in the quasi-flat part of the light curve, it should be 
possible to estimate the density of the cloud, or at least 
an upper limit on it, to within a factor of 2. The second 
(later) measurement is desirable because it enables us to 
roughly constrain the shape factor f{t). To determine the 
collimation angle it is necessary to have a better esti- 
mate of ti, which requires at least three flux measurements 
at different times to be carried out, and also to know the 
radius of the cloud R. One can estimate the latter quan- 
tity, which is of interest in its own right, by measuring 
the absorption column density toward the GRB or its af- 
terglow: R = Nn/n. It is of crucial importance that the 
dependence of Oq on all parameters is fairly weak - see 
equation (pi]). 

The actual environment of a GRB may differ signif- 
icantly from our simplified model. However, we expect 
that the shape of the light curve of the scattered emis- 
sion will not be affected dramatically by this. To illustrate 
this point we have plotted in Fig. ^ three different com- 
puted light curves corresponding to the case of 9q = 25° 
and R = 10 pc, one of which is described by our ana- 
lytic solution. We see that these curves are very similar 
despite the significant differences in the setup. This has a 
straightforward explanation: the signal detected at a given 
time results from integration over scattering sites, satisfy- 
ing the condition 9i < 9 < Oa, that are located at different 
distances and in different directions from the center of the 
cloud. Therefore, the scattered flux is proportional to some 
weighted average, (n) , of the gas density over the scatter- 
ing surface, which only weakly responds to changes in the 
geometry and distribution of gas. 

Although we have so far only considered the scat- 
tering process, the GRB radiation can also be photoab- 
sorbed by the neutral (or weakly ionized) matter, giving 
rise to fluorescence emission, mainly in the Fe Ka line 
(Ghisellini et al. 1999; Bottcher et al. 1999). Given the 
dependence of the photoabsorption cross section on en- 
ergy, cr w 3.5 X 10-2"cm2(7.11kcV/£;)2•^ and the fluo- 
rescence yield of 0.34 (e.g. Vainshtein & Sunyaev 1981), 
assuming that the spectrum of primary X-ray emission is 



dNphoton/dE oc i?-^ with 7 close to 1 (as typically ob- 
served, see e.g. Amati 2002), and assuming solar abun- 
dance of iron, we find that fluoresence radiation con- 
tributes a fraction ~ 0.2A~^{1 + cos^ 9)~^ to the total 
detectable flux at 0.3-5 keV (which corresponds to the 
rest-frame 0.6-10 keV at z = 1), i.e. ^ 7% for the most 
interesting angles 9^0 and a molecular fraction A ^ 1.5. 
This estimate is valid so far as the X-ray spectrum is not 
heavily absorbed, i.e. for A^h ^ a few xlO^^ cm-^, and is 
not strongly dependent on 7 and z. We may therefore con- 
clude that the fluorescence emission, although interesting, 
will typically be an order of magnitude fainter than the 
scattered component. 

3. Detectability of the effect 

A key question is of course the detectability of the X- 
ray echo. Equation (^) indicates that the effect should 
be just within the reach of the currently flying Chandra 
and XMM-Newton observatories for bursts with S'x ~ 
10-5 erg cm-2 occuring in molecular clouds with n ^ 
10"^ cm-^. How realistic are these values for both param- 
eters? 

First, we know that GRBs having Sx ^ 10-5 ergcm-^ 
do occur from time to time. Note that it is the X-ray 
fluence, rather than the total fluence, which is important 
for us. To make an estimate we have taken information 
(Amati et al. 2002) on the fluxes and spectra at > 2 keV 
of twelve bright GRBs observed with the BeppoSAX satel- 
lite, then extrapolated the spectral fits to lower energies 
and calculated the fluences in the 0.3-5 keV band, where 
both Chandra and XMM-Newton have good sensitivity. 
We find two bursts among this sample with a sufficiently 
high X-ray fluence: GRB 990712 (S^x = 5 x 10-^ erg cm-^) 
and GRB 010222 (^x = 7 x 10-^ erg cm-^). Other 
gamma-ray burst detectors with sensitivity in the X-ray 
band such as the one on Ginga (Yoshida et al. 1989), 
GRANAT/WATCH (Sazonov et al. 1998) and recently 
HETE/FREGAT (Barraud et al. 2002) have also detected 
a few GRBs with Sx ^ 10-5 erg cm-^. 

Fluence values quoted in the literature usually pertain 
to the prompt GRB emission. However, there are indica- 
tions (Burenin et al. 1999; Giblin et al. 1999; Tkachenko 
et al. 2000) that a comparable X-ray fluence may be con- 
tained in early (~ 10^ s) GRB afterglows, during which 
the energy spectrum is much softer than during the burst 
proper. Without attracting spectral information it cannot 
be possible to separate the contributions of the prompt 
and early afterglow emission to the X-ray echo at early 
times {t <C ii), as the scattered X-ray flux will be sim- 
ply proportional to the total of the X-ray fluences of the 
burst proper and its early afterglow. However, there are 
two possibilities for the light curve of the X-ray echo. In 
the case where the early afterglow emission is beamed ex- 
actly as the prompt radiation, the conclusions of all 
remain true, except that in equation (^) one should un- 
derstand Sx as the combined X-ray fluence of the GRB 
and the early afterglow. If, on the other hand, the early af- 
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terglow is characterized by its own beaming angle ^ Oo 
the light curve of the X-ray echo will be a sum of two 
similarly shaped light curves (like those shown in Fig. 
one corresponding to the other to with their 

relative weights being proportional to the fluences of the 
prompt and early afterglow X-ray emission. Therefore, it 
should be in principle possible to estimate from the scat- 
tered light curve the angle 6'a in addition to Oq. 

As regards the amplitude of the effect, we should also 
mention that Galactic absorption will lead to a reduction 
of the GRB fluence and similarly of the scattered flux in 
the 0.3-5 keV band. This, however, will only affect the flux 
near the lower boundary of the quoted spectral range, so 
the net effect is expected to be small 10%). 

Let us next consider the gas density. This is of course 
a quantity which the proposed effect enables to constrain. 
We may, however, speculate a little given the information 
available now. As mentioned in §|l|. X-ray spectral mea- 
surements reveal substantial optical depths to photoab- 
sorption, Nn ^ 10^^ cm^^, in the directions of GRBs. 
Such column densities are similar to those of giant molec- 
ular clouds in the Milky Way and therefore imply number 
densities 71 - 3x10^ cmr^{NulW^ cm-2)(io pc/i?). This 
is about what is needed for the scattered X-ray emission 
from the brightest GRBs to be detectable. There is ongo- 
ing debate (e.g. Galama et al. 2001) as to why measured 
optical extinctions tend to be small even when the corre- 
sponding A^H is large. A possible explanation is that the 
UV and X-ray radiation from the GRB and its afterglow 
evaporates dust out to a few times 10 pc (Waxman & 
Draine 2000; Fruchter et al. 2001). We should also men- 
tion the issue of "dark bursts" , i.e. GRBs with undetected 
optical afterglows. It is possible that such bursts occur in 
very dense clouds so that their optical afterglow emission 
is extinguished (e.g. Reichart & Yost 2001; Ramirez-Ruiz 
et al. 2002). 

It is interesting to consider here again (see §||) the 
case of GRB 000926, whose afterglow was observed by 
BeppoSAX and Chandra. From analysis of X-ray and op- 
tical data, Piro et al. (2001) have inferred that the jet had 
an opening angle of 0o ~ 25° and expanded into a dense 
medium with n = 4 x 10*^ cm^'^ (see, however, Harrison 
et al. 2001; Panaitescu & Kumar 2002). At the time of 
the latest Chandra observation [t ~ 10^ s), the afterglow 
emission was still easily detectable, with the flux at 0.2- 
5 keV being 8 x 10" gj.g ^^-2 g-i^ Qj^g 000926 was 
discovered by the Interplanetory Network (Hurley et al. 
2000), and its 25-100 keV fluence was 6.2 x 10"^ erg cm'^ 
(Piro et al. 2001). Assuming different values for the slope 
7 of the X-ray part of the GRB spectrum we may es- 
timate that the GRB 000926 fluence at 0.3-5 keV was 
4 X 10""^, 2 X 10^^ erg cm^^ and 1.2 x 10~^ erg cm"^ 
for 7 = —1, —1.5 and —2, respectively. We can then esti- 
mate from equation (|^) the scattered X-ray flux during the 
last Chandra observation (assuming that t <^ ti so that 
f{t) - a few): Fx - 10"i^(l + z)-\ 4x 10"i''^(l + z)"i and 
2 X 10~^'*(1 -I- z)~^ erg cm~^ s~^ for the 7 values given 
above. Therefore, if the density of the medium around 



GRB 000926 is indeed as high as implied by the analysis 
of Piro et al. (2001), then the scattered flux already could 
have constituted a significant fraction of the observed X- 
ray flux two weeks after the burst. In this case, however, 
the size of the scattering cloud must not be too large 
{R <^ 10 pc). Otherwise, it would be difficult to explain the 
marginally detected column density A'h ~ 4 x 10^^ cm^^ 
(Piro ct al. 2001), even taking into account the pho- 
toionization effect of the GRB on the ambient medium 
(Bottcher et al. 1999). 

4. Conclusions 

We have described a new observational method that en- 
ables us to determine the typical ambient medium density 
of bright GRBs as well as their degree of collimation in an 
almost model-independent way. This method can hope- 
fully be employed with the existing Chandra and XMM- 
Newton X-ray telescopes. 

A suitable observational strategy would be to sched- 
ule an observation of the site of a very bright (5'x ~ 
10~^ erg cm~^) GRB a few months after the burst. Should 
some X-ray flux be detected during this observation, an 
additional one or two observations should be carried out 
several months or years later so that the light curve of 
the scattered emission could be roughly reconstructed and 
the GRB opening angle could be estimated [from equa- 
tion (|ll|)]. If the first observation fails to yield a signifi- 
cant flux, an interesting upper limit, ^ 10'^ cm~'^, on the 
density of the medium surrounding the GRB on parsec 
scales can be derived [from equation (|l^)]. It will take a 
10^ s observation with Chandra or XMM-Newton to de- 
tect a scattered flux of a few xl0~^^ erg cm~^ s""'^ at 
0.3-5 keV. We emphasize that it is sufficient to collect a 
total of several photons from the GRB direction, as no 
detailed spectral information is needed. We also note that 
a flux of 5 X 10~^^ erg cm~^ s~^ corresponds to an equiv- 
alent isotropic luminosity of order lO''^ erg s~^ for a burst 
at z ~ 1. Therefore, the scattered GRB emission should 
outshine the X-ray emission of any non-active host galaxy. 

Since GRBs with S'x ~ 10^^ erg cm~^ are very rare 
events (at best a few per year), constant monitoring of the 
whole sky with an instrument sensitive to X-rays is crucial. 
After the end of the CGRO mission, such a capability will 
be provided by the Swift satellite^ 
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